Probing anisotropic superfluidity of rashbons in atomic Fermi gases 
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Motivated by the prospect of realizing a Fermi gas of ^''K atoms with a synthetic non-Abelian 
gauge field, we investigate theoretically a strongly interacting Fermi gas in the presence of a Rashba 
spin-orbit coupling. As the two-fold spin degeneracy is lifted by spin-orbit interaction, bound pairs 
with mixed singlet and triplet pairings (referred to as rashbons) emerge, leading to an anisotropic 
superfiuid. We show that this anisotropic superfiuidity can be probed via measuring the momentum 
distribution and single-particle spectral function in a trapped atomic "'"K cloud near a Feshbach 
resonance. 
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Owing to the unprecedented experimental control- 
labihty, ultracold atoms have been proven to be an 
ideal tabel-top system to study some long-sought, chal- 
lenging many-body problems. A well-known example 
is the non-perturbative problem of crossover from a 
Bose-Einstein condensation (BEG) to a Bardeen-Cooper- 
SchriefFer (BCS) superfluidity in an ultracold atomic 
Fermi gases Here we study a strongly interacting 
Fermi gas in the presence of a synthetic non-Abelian 
gauge field, as motivated by the recent demonstration 
of such field in bosonic ^^Rb atoms and the prospect 
of its realization in fermionic "^^K atoms We focus on 
the Rashba spin-orbit (SO) interaction and explore its 
impact on the unitary Fermi gas. 

A weak non-Abelian gauge field such as Rashba SO 
interaction is well understood. It was shown in 2001 by 
Gor'kov and Rashba that a condensed matter system 
of superconducting 2D metals with weak SO coupling fe- 
tures a mixed spin singlet-triplet pairing field, and its 
spin magnetic susceptibility can be dramatically affected 
by the SO interaction. By applying an additional large 
Zeeman magnetic field, it was proposed by Zhang et al. 
and Sato et al. [g that a topological phase with 
gapless edge states and non-Abelian Majorana fermionic 
quasiparticles may form. More recently, Vyasanakere and 
Shenoy identified an interesting bound state by solving 
the two-body problem Q, referred to as rashbons. By 
increasing the strength of SO coupling, a BCS superfiuid 
can therefore evolve into a BEG of rashbons Q. 

In this Letter, we investigate the salient features of 
rashbons by using a functional path-integral functional 
method. We identify clearly rashbons from the gaus- 
sian fluctuations of pairing field and show they possess 
anistropic effective mass. We demonstrate at the saddle- 
point level that the condensation of rashbons gives rise to 
both singlet and triplet pairings. As a result, the momen- 
tum distribution and single-particle spectral function be- 
come highly anisotropic. By performing calculations for 
a realistic system of trapped '^^K fermions at a Feshbach 



resonance, we present observable experimental signatures 
for visualizing the anisotropic superfluidity of rashbons. 

The model — Let us start by formulating the BEG- 
BCS crossover with a Rashba SO coupling T-Lso = 
^{kydx — kx(Ty), whose Hamiltonian is given by, 

H= j dv {v>+ + Uso] + Uo^+^+^^^^] , (1) 

where = h^k^ / {Im) — /i. and V'(r) = [iI]^{y),iP:^{y)] 
denotes collectively the fermionic field operators. 
The contact s-wave interaction (C/q < 0) occurs 
between un-like spins. We use the functional 
path integral method Q and consider the action 
Z = i/i] exp {— 5 [?/; (r, r) (r, r)] }, where 

S[ij,i,] = jyrij dvj:^j;, (r) dri^^r)] + H{i:,i;), 
/? = 1/(/cbT), and Hlipji/j) is obtained by replacing 
the field operators ■ip'^ and ip with the Grassmann 
variables tp and ip, respectively. The interaction term 
can be decoupled by using the standard Hubbard- 
Stratonovich transformation with the introduction 
of a fluctuating pairing field A(r,T). It is con- 
venient to use the 4-dimensional Nambu spinor 
(f> (r,T) = [V't, V'^iV'ti '04,] a-iid rewrite the action as, Z = 
/ P[<I>, $;A, A] exp{- /^^ dr J dv [i$ {-g-^) $ - 
/^Sk^k}, where Q is the single-particle Green function. 

By integrating out the fermionic field, the effective ac- 
tion isZ = J V[A, A] exp {-5'e// [A, A] }, where Seff = 

/„^rfr/dr{-^(lii^}-iTrln[-g-i]+/?Ekek. To 
proceed, we restrict ourselves to the gaussian fluctua- 
tion and expand A(r,r) = Aq + 5A(r,T). The effec- 
tive action is then approximated by 5'e// = Sq + AS, 

where the saddle-point action 5o = f^^drfdr^^—^ 

^Trln [—0(7^] + /3 X^k and, in the momentum space, 
the fluctuation action takes the form [k = (k, iwm) 

and q = (q,w„)]: AS = 
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Rashbon — Let us consider first the normal state with 
Aq = 0, in which case the Green function reduces to its 
non-interacting form as Go{k) = Diag{go(fc)7 ~3o(~fc)} 
with gaik) = [iuj^ - - X{kya^ - k^cTy)]'''-, lead- 
ing to two helicity branches in the single-particle spec- 
trum, i?k,a = + a\k±, where k± = (k^ + kyY^^ and 
a = ±1. The fluctuation action is given by AS = 
[-r^i (q)] 5A{q)SA{q), where T^'^iq) is the inverse 
vertex function which, at q = 0, takes the form 



E 

a=± 



1/2 -/(gk,^ 



1 

2^ 



where f (x) = Xjie^l^'^'^ -I- 1) is the Fermi distribution 
function and we have renormalized the bare interaction 
C/q by the s-wave scattering length, 1/f/o ~ m/(47rfi^as) — 
'YliY l/(2ek), with V being the quantization volume. 
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Figure 1: (Color online) (a) Raslibons as evidenced by the 
two-body phase shift of F^^ (0,a;) at three different scatter- 
ing lengths. The arrows indicate the position of binding en- 
ergy, (b) Effective mass of rashbons [7 = Mx/(2m)] in the 
strong SO hmit. The inset shows the bound state energy as 
a function of the scattering length. 

The vertex function is simply the Green function of 
the fermion pair. A bound state can therefore be exam- 
ined clearly by calculating the phase shift i5(q, w) = 
— Imln[— F~^(q, w„ lo + iO+)]. For a true boson, the 
phase shift is given by (S_B(q, lo) — TTQ{uj~e^+fiB), where 
Eq and hb are the bosonic dispersion and chemical po- 
tential, respectively, and ld{x) is the step function. In 
Fig. [Ija), we plot the two-body part of the phase shift 
at q = 0, obtained by discarding Fermi functions. The 
phase shift jumps from to tt at a critical frequency, re- 
sembling that of a true boson. This is exactly the demon- 
stration of a bound state, or rashbon. By recalling that 



the bosonic chemical potential is given by /zs = 2/x — -Bs, 
where Eb is the bound state energy, the critical frequency 
{uj + 2/.j)c at q = gives exactly Eb- Using the fact that 
the critical frequency corresponds to the position where 
Rc[F^^] changes sign, we have 
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The inset in Fig. [Ijb) shows the bound state energy as 
a function of the SO coupling strength. At the unitar- 
ity limit, the bound state energy is universally given by 
EB{a, = ±00) « -1. 439229mA V^^- The size of rash- 
bons a is therefore at the order of h?/ (to A). Rashbons are 
well defined once a <^ kp^ or Xkp 3> ef- Thus, we antic- 
ipate that the system will cross over to a gas of rashbons 
at Xkp/iF ~ 1- 

In the limit of a large SO coupling, the well- 
defined rashbons should have a bosonic dispersion = 
h^q\/{2Mj^) + h^ql/{2M^) and weakly interact with each 
other repulsively. Because of the anisotropic fermionic 
dispersion -Bk.i = ■Ck ± Afcj^, the effective mass of rash- 
bons becomes anisotropic. While ~ 2m is not af- 
fected by the Rashba coupling, AI± may get strongly 
renormalized. As the jump of the phase shift at nonzero 
q which occurs at (w -|- 2fj,)c — Cq , we can numerically 
determine M±. Fig. [TJb) reports 7 = AI±/{2m). At uni- 
tarity, we find 7 ~ 1.2. When the system becomes an 
ensemble of weakly interacting rashbons, the heavy mass 
M± causes a decrease in the condensation temperature 
so that Tbec = 7~^^^rBEC,o, where Tbec.o - 0.218Tf is 
the BEG temperature without the SO coupling. 

Condensation of rashbons — Let us now turn to the 
condensed phase characterized by a nonzero order pa- 
rameter Aq 7^ 0. At the mean-field saddle-point level, 
the single-particle Green function takes the form. 



Go'- 



, - a - Us 
-iAo&y 



iAoo-y 



(4) 



The inversion of the above matrix can be worked 
out explicitly, leading to two single-particle Bogoli- 
ubov dispersions whose degeneracy is lifted by the 

SO interaction, -Ek,± = \J (Ck ± Afc^)^ + Aq, and the 
normal and anomalous Green functions from which 
we can immediately obtain the momentum distribu- 
tion n (k) = 1 - J2a [1/2 - / {Ek,a)] 7k, a and the 
single-particle spectral function ^-[-(k, cj) = ^^^(k, w) = 

Ea [(1 + 7k,a) S{lU~ £;k,a) + (1 " 7k,a) <5 (w + E^^^c.)] /4, 

where 7k, ± = (Ck i ^k±^) /-Ek,±- The chemical potential 
and the order parameter are to be determined by the 



number and the gap equations. 



X]k"-(k) and 



Ao = -t/oAoEa[l/2 - /(^^k,a)]/(2£;k,a), respectively. 
Fig. [2] displays the chemical potential and order pa- 
rameter as functions of the SO coupling strength for 
a unitary Fermi gas. The increase of the SO strength 
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leads to a deeper bound state. In analogy with the 
BEC-BCS crossover, the order parameter and critical 
transition temperature are greatly enhanced once the 
rashbons are well defined at Xkp ^ ep. In the large SO 
coupling limit, the chemical potential can be written 
as fjL ~ Eb l'^ + yi-B I'^i where the chemical potential of 
rashbons //s is positive due to the repulsion between 
rashbons and decreases with increasing coupling as 
shown in the inset of Fig. EJ^c). By assuming an 
s-wave repulsion with scattering length as, where 
\ji,B — {n/'2)Ai:h^aB / M , we estimate within mean-field 
that in the unitarity limit, as — 3h'^/{mX), comparable 
to the size of rashbons. A more accurate description of 
the interaction between two rashbons can be obtained 
by solving a four-fermion problem. 




Figure 2: (a) Momentum distribution and (b) single-particle 
spectral function for 6 — n/2 at Xk^/ep = 2. Here 6 is the 
angle between k and the z-axis. The width of the curves in 
(c) represents the weight factor (1 ±7k,±)/4 for each of the 
four Bogoliubov excitations, (c) Mean-field order parameter 
as a function of the SO coupling for a homogeneous unitary 
Fermi gas at zero temperature. The inset shows the chemical 
potential and the half of bound state energy, in units of ef- 

Figure [H^a) and (b) illustrate the momentum distri- 
bution and the single-particle spectral function, respec- 
tively. These quantities exhibit anisotropic distribution 
in momentum space due to the SO coupling and can be 
readily measured in experiment. 

Another important consequence of the SO cou- 
pling is that the pairing field contains both a sin- 
glet and a triplet component Q, Si- For the sys- 
tem under study, it is straightforward to show that 
the triplet and singlet pairing fields are given by 
iM-n) = -*Aoe-''^>'Ea«[l/2 - /(i^k.a)]/(2£;k,a) 
and (V'ktV'-ki) = AoX;Q[l/2-/(-Bk,a)]/(2£'k,a), respec- 



tively, where 



(fca; — iky)/kj_. The magnitude of 



these are shown in Fig. |31 The weight of the triplet com- 
ponent increases as the SO coupling strength increases 
and becomes equal to that of the singlet component in 
the strong SO coupling limit. 

Probing the anisotropic superfluidity of rashbons — 
One leading candidate to observe superfiuid rashbons is 
a trapped Fermi gas of atoms near a broad Fesh- 
bach resonance, where an applicable scheme to generate 
the Rashba SO coupling was recently proposed |3|. Pre- 
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Figure 3: (Color onhne) Linear contour plot for the triple pair- 
ing correlation | (V'kt'/'-kt) I between like spins (upper panel) 
and the singlet pairing correlation | (i/'ktV'-ki) | between un- 
like spins (lower panel) for a homogeneous unitary Fermi gas 
at zero temperature. The SO coupling strength in (a) and (b) 
is Xkp/ep = 1, and in (c) and (d) is Xkp/ep — 2. 
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Figure 4: (Color online) Density distributions of a trapped 
unitary Fermi gas at zero T = (a) and finite temperature 
T — O.STc.o (b) and at different SO coupling strengths. Tc.o 
is the critical temperature in the absence of the SO coupling. 
The inset in (a) shows log-scale contour plot of the momen- 
tum distribution and the inset in (b) illustrates the profile of 
the order parameter at finite temperature. Here the Fermi en- 
ergy is given by Ep = (SNy^^hujo, Fermi wave number kp = 
(24iV)^''^a,7j, Thomas-Fermi radius ttp = (24iV)^''^aho, and 
the non-interacting peak density upp = {24N)^^^ /(3Tv''^)a 
where aho ~ y/h/{muj) is the characteristic length of the har 
monic oscillator. 
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vious experiments have demonstrated the measurement 
of momentum distribution and single-particle spectral 
function is '^'^K in the absence of the SO coupling 
We perform the mean-field calculation in a 3D spheri- 
cal harmonic trap Vrir) = mujQr'^/2, by using the lo- 
cal density approximation (LDA). In LDA, the gas was 
divided into small cells with a local chemical potential 
/i(r) = /i — Vrir). The local density n(r), momen- 
tum distribution n(k; r)/(27r)^, occupied spectral func- 
tion ^(k, uj; r)f {(jj) fe^/ (27r^) are then integrated over the 
whole space to obtain the total contribution. The chem- 
ical potential at the trap center fi can be determined 
by the number equation J drn{r) = N, where N is the 
total number of fermions. We show, in Fig. |4l the den- 
sity profile of a trapped unitary Fermi gas at different 
SO coupling strengths and temperatures. As anticipated, 
with the increase of the SO coupling the cloud shrinks. 
The anisotropic momentum distribution at large SO cou- 
pling, which can be measured using the time-of-flight 
technique, can still be clearly seen as illustrated in the 
inset of Fig. [H^a) . 
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Figure 5: (Color online) Log-scale contour plot of the single- 
particle spectral function for a trapped unitary Fermi gas at 
T = O.STco- The SO coupling strength in (a) and (b) is 
Xkp/EF — 1, and in (c) and (d) is Xkr/Ep = 2. 

Figure [S] presents the occupied spectral function at 
T ~ O.STc.o, where is the critical temperature with- 
out the SO interaction. The distinct behavior for the 
spectral function at 6* = (along the z-axis) and tt/2 
(in the transverse plane) can be probed by the rece ntly 
demonstrated momentum-resolved rf-spectroscopy 
To make better comparison with the experiment, we have 
included in the calculation an energy resolution of 0.2Ep, 
as presented in the JILA rf- measurement. Our calcula- 
tions therefore clearly demonstrate that the anisotropic 
nature of the rashbon superfluid will not be smeared out 
by averaging over the trapped cloud. 

Summary — We have shown that rashbons induced by 
a strong Rashba spin-orbit coupling differ significantly 
from the BEC-BCS crossover molecules studied over the 



past few years. Rashbons have anisotrpic effective mass 
and condense into a mixed singlet-triplet pairing state. 
They lead to a strong anisotropy in the momentum dis- 
tribution and the single-particle spectral function. We 
have proposed that these distinct behaviors can be read- 
ily probed in a strongly interacting Fermi gas of "^^K 
atoms in a synthetic non-Abelian gauge field. 

More interesting features of rashbons may be discov- 
ered thanks to the unprecedented controllability in ul- 
tracold atoms. Rashbons in 2D may be utilized to cre- 
ate Majorana fermions . Rashbons under a mag- 
netic field may exhibit inhomogeneous superfluidity with 
mixed singlet-triplet pairings and their properties may be 
greatly affected by the field. Therefore, the exploration 
of strong correlation effects of rashbons is a new excit- 
ing many-body problem. In the current work, we have 
adopted a mean-field approach. The mean-field calcula- 
tion in the condensed phase can be improved by incorpo- 
rating gaussian pair fluctuations (l3| . 
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Note added: After we completed our work, we be- 
came aware of the work reported in Ref. |14l | which treats 
a similar system. Our results agree with each other where 
they overlap. 
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